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We show that the plasmon spectrum of an ordinary two-dimensional electron gas (2DEG) hosted 
in a GaAs heterostructure is significantly modified when a graphene sheet is placed on the surface of 
the semiconductor in close proximity to the 2DEG. Long-range Coulomb interactions between mas- 
sive electrons and massless Dirac fermions lead to a new set of optical and acoustic intra-subband 
plasmons. Here we compute the dispersion of these coupled modes within the Random Phase Ap- 
proximation, providing analytical expressions in the long-wavelength limit that shed light on their 
dependence on the Dirac velocity and Dirac-fermion density. We also evaluate the resistivity in 
a Coulomb-drag transport setup. These Dirac/Schroedinger hybrid electron systems are experi- 
mentally feasible and open new research opportunities for fundamental studies of electron-electron 
interaction effects in two spatial dimensions. 
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I. INTRODUCTION 

Plasmons are ubiquitous high-frequency collective den- 
sity oscillations of an electron liquid, which occur in met- 
als and semiconductors^. Their importance across differ- 
ent fields of basic and applied physics is by now well 
established. They play, for example, a key role in plas- 
monics^ and in the photodetection of far-infrared radia- 
tion based on field-effect transistors^. 

Plasmons in doped graphene sheets (hereafter dubbed 
"Dirac plasmons") have been intensively investigated 
by electron energy-loss spectroscopy (EELSpl EELS, 
however, does not have the sufficient energy and wave- 
number resolution to address e.g. the impact of bro- 
ken Galilean invariance on Dirac plasmons^. Moreover, 
its application is limited to the collective modes in the 
charge channel with no information on the spin degrees 
of freedom, since the method probes only the "loss" 
function, i.e. the imaginary part of the inverse di- 
electric function e(q,ui). Angle-resolved photoemission 
spectroscopy^ and scanning tunneling spectroscopy^ also 
give precious information on Dirac plasmons, albeit in 
a slightly-indirect way. Finally, Dirac plasmons have 
been probed by engineering directly their coupling to far- 
infrared light in a number of intriguing ways^. On the 
contrary, the application of inelastic light (Raman) scat- 
tering to probe plasmons in graphene, has been, so far, 
limited to the investigation of magneto-plasmons in high 
magnetic fields^. 

The situation is very different for the case of plas- 
mons in 2DEGs realized, for example, in GaAs/AlGaAs 
modulation-doped semiconductor heterostructures. In- 
deed, these excitations have been successfully studied, 
both in zero magnetic field and in the qu antum Hall 
regime, by resonant inelastic light scatterin g* 11 ! 12 . This 
method has a much better energy resolution than EELS 
(typically below 0.1 meV) and can be straightforwardly 
generalized to probe spin excitations. The cross section 
for scattering of photons by the electronic subsystem is, 
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FIG. 1: (Color online) Schematics of the double-layer mass- 
less Dirac/Schroedinger hybrid electron system studied in 
this work. A graphene sheet is deposited on the surface of 
a semiconductor, underneath which a GaAs quantum well 
(at a distance d from the surface) hosts a high-mobility two- 
dimensional electron gas (2DEG). Carriers in the two layers 
are induced by conventional gating techniques (the + signs 
below the quantum well indicate the doping layer). In the 
presence of Coulomb coupling between the two subsystems, 
optical and acoustic hybrid collective modes emerge, which 
can be probed by resonant inelastic light scattering. 



however, tiny. An ingenious resonance condition can be 
used in a light scattering experiment to enhance this cross 
section dramatically^. The frequency of the incoming or 
scattered photon can indeed be tuned close to the semi- 
conductor band gap {e.g. « 1.5 eV for GaAs at low 
temperature) making the energy denominators appear- 
ing in second-order perturbation theory^- 3 small and thus 
yielding a much higher sensitivity to collective electronic 
modes. 

Graphene^ is a gapless material though, and thus 
presents immediately an obstacle for the application of 
resonant inelastic light scattering. The van Hove singu- 
larity at the M point of graphene's Brillouin zone can 
be used to achieve resonance^, although the band gap 
at the M point is in the ultraviolet region of the electro- 
magnetic spectrum (« 4 eV) where tunable lasers having 
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the required small linewidth for resonant inelastic light 
scattering do not exist. 

In this Article we propose a hybrid double-layer system 
composed by a doped graphene sheet that is Coulomb 
coupled to an ordinary 2DEG in a GaAs/AlGaAs semi- 
conductor heterostructure, as schematically illustrated 
in Fig. [T] We demonstrate that precious information 
on Dirac plasmons is encoded in the dispersion of the 
new collective in-phase and out-of-phase plasmon exci- 
tations of this massless Dirac/Schrocdinger hybrid sys- 
tem (MDSHS). These collective modes can in principle 
be probed via inelastic light scattering by exploiting, as 
explained above, the resonance condition offered by the 
GaAs band gap. In passing, we mention that a simi- 
lar resonant enhancement mechanism has been recently 
observed in the case of the Raman phonon lines in epi- 
taxial graphene on copper substrates with incident pho- 
ton energies overlapping the copper substrate photolumi- 
nescenceiS. We also study transport in a MDSHS in a 
Coulomb-drag setup^, demonstrating that, in the low- 
temperature Fermi-liquid regime and at strong coupling, 
the drag resistivity displays an intriguing dependence on 
carrier densities and inter-layer separation, which is dif- 
ferent from that knowrf^ for coupled 2DEGs in semicon- 
ductor double quantum wells. 

Our manuscript is organized as follows. In Sect. [H] we 
introduce the model Hamiltonian and some basic defi- 
nitions. In Sect. |IIl| we present analytical and numerical 
results, obtained by means of the Random Phase Approx- 
imation, for the collective plasmons modes of MDSHSs. 
Our main analytical results on the Coulomb drag resis- 
tivity are summarized and discussed in Sect. |IV] Finally, 
in Sect. [V] we present a short summary of our main results 
accompanied by a discussion of the potential of MDSHSs. 



II. MODEL HAMILTONIAN AND IMPORTANT 
DEFINITIONS 

The system depicted in Fig. [I] can be modeled by the 
following Hamiltonian: 
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The terms in the first and second lines of the previous 
equation describe two 2D electron systems with massless 
Dirac and parabolic Schroedinger bands, respectively. 
The last term describes Coulomb interactions between 
the two subsystems. 

In the first line of Eq. ([lj, «rj is the Dirac velocity, 
er = (a x ,a y ) is a 2D vector of Pauli matrices, S is the 



sample area, and %j)' h (ipk.a) creates (destroys) an elec- 
tron with momentum k and sublattice-pseudospin index 
a = A, B in graphene. The Dirac density operator is de- 
fined by ffi = J2k,a ^k-q,a^,a- In the second line, m h 
is the bare band mass of the 2DEG (mb ~ 0.067 m c in 
GaAs, m e being the electron mass in vacuum), $ k (</>&) 
creates (destroys) an electron with momentum k in the 

2DEG, and p q h) = £ fc 4>k- q 4>k is the usual 2DEG density 
operator. 

In writing Eq. we have intentionally hidden "flavor" 
indices labeling spin and valley degrees-of-freedom in the 
graphene part of the Hamiltonian and spin degrees-of- 
freedom in the 2DEG part of the Hamiltonian. These 
flavors will just appear as degeneracy factors in the lin- 
ear response functions xt(?,w) and Xh{q,u), as we will 
see below. We have also completely neglected inter-layer 
tunneling, which is strongly suppressed by the AlGaAs 
barrier between graphene and the quantum well. 

The bare intra- and inter-layer Coulomb interactions, 
Vij(q) (with i,j £ {t,b}), are influenced by the lay- 
ered dielectric environment, modeled by two dielectric 
constants - ei and e 2 in Fig. [T] A simple electro- 
static calculatio n 1 18 ! 19 ! implies that the Coulomb inter- 
action in the graphene sheet (top layer) is given by 
Vu(q) = 47re 2 /[(?(ei + e 2 )], The Coulomb interaction in 
the 2DEG (bottom layer) is instead 
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where D(q) = 2e 2 (ei + e 2 )e 9 
interaction is given by 

V th {q) = V ht (q) 



Finally, the inter-layer 
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62 = 6 limit we recover 

2 , 
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Note that in the "uniform" t\ 
the familiar expressions Vu(q) = Vbb(q) 
V th (q) = V ht (q) = V tt (q) exp(-qd). 

For the following analysis we introduce the electron 
density in the top (bottom) layer n t (rib) and the corre- 
sponding Fermi wave numbers /cf,i = \J 47rrij ~/Ni , where 
Ni is a degeneracy factor (N t — 4 in the graphene 
layer and N h = 2 in the 2DEG layer). We also in- 
troduce the Fermi energies £F,t = fivukF,t ano ^ £ F,b = 
h 2 kp b /(2mb) in the top and bottom layers, respec- 
tively, a oc = e 2 /(hv D ) fa 2.2, and r s = (7m b a|)~ 1/2 , 
ae = h 2 /(mbe 2 ) being the Bohr radius calculated with 
the semiconductor band mass nib (note the absence of 
any dielectric constant in the definition of ae). 



III. COLLECTIVE OPTICAL AND ACOUSTIC 
PLASMON MODES 



The collective modes of the system described by the 
Hamiltonian (TIT) can be determined^ by locating the poles 
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of the linear-response function x(Qi u )- Within the Ran- 
dom Phase Approximation (RPA) these functions satisfy 
the following matrix equation^ 



- 1 (q,tu)=Xo 1 (q^)-V(q) , 



(4) 



where Xo(<Zi w) is a 2 x 2 diagonal matrix whose elements 
Xtiljto) and Xb(<?j w ) arc the well-known non-interacting 
(Lindhard) response functions of each layer at arbitrary 
doping m. The mathematical and physical properties 
of Xt(<7> w ) are discussed in Ref. [201 while expressions 
for 5Re [xb(9, w )] and [xbfew)] can be found in 
Ref. [TJ The off-diagonal (diagonal) elements of the ma- 
trix V = {Vij}i t j—t,b represent inter-layer (intra-layer) 
Coulomb interactions. 

A straightforward inversion of Eq. ^ yields the fol- 
lowing condition for collective modes: 



e(q,u) = [l-Vtt(q)Xt(q,w)][l- 
- V? h (q)xt{q,u)Xb(q,u) 



Vbb(q)xb(q,w)] 

= 0. (5) 



Zeroes of e(q, uj) occur above the intra-band particle-hole 
continuum where Xi is real, positive, and a decreasing 
function of frequency. Eq. ^ admits two solutions^, 
a higher frequency solution at ui op (q — > 0) oc Jq which 
corresponds to in-phase oscillations of the densities in 
the two layers (optical plasmon), and a lower frequency 
solution at w ac (g — > 0) oc q which corresponds to out-of- 
phase oscillations (acoustic plasmon). 

A summary of our main results for w op (g), obtained 
from the numerical solution of Eq. ([5]), is reported in 
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In these figures we plot the quantity A(q) 
huj p i(q), which physically represents the energy 
of the MDSHS optical plasmon measured from that of 
the plasmon of an isolated 2DEG, hui p \(q). We clearly 
see that A(q) is positive and « 10 meV at q ~ 10 5 cm -1 . 

To leading order in q in the long-wavelength q — > 
limit, the frequency of the MDSHS optical plasmon mode 
can be found analytically. The result is 
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where we have introduced e = (ei + e 2 )/2. We have 
checked that Eq. |6| is in perfect agreement with the 
numerical results displayed in Figs. 2]|3 The first term 
on the r.h.s. of Eq. can be easily recognized^ to be 
the square of the RPA plasmon frequency of the elec- 
tron gas in an isolated graphene sheet separating two 
media with dielectric constants t\ and e 2 . The sec- 
ond term is the well-known RPA plasmon frequency of 
a 2D parabolic-band electron gaffi A measurement of 
u! p(q) m a MDSHS at sufficiently small q [more pre- 
cisely, for q <C min(A:F j t, fcp,b)] thus allows to access di- 
rectly the Dirac velocity «d and the Dirac-fermion den- 
sity n t cx kp t . 

In Fig. [4] we show illustrative results for the acous- 
tic plasmon w ac (g) as obtained from the numerical so- 
lution of Eq. ([5]) for different values of the 2DEG car- 
rier densities n^. The dashed line in this figure repre- 
sents the upper bound u> = vj^q of the Dirac-fermion 
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FIG. 2: (Color online) Optical plasmon dispersion in a mass- 
less Dirac/Schroedinger hybrid electron system. In this fig- 
ure we plot the quantity A(q) = hjj op (q) — hw p i(q) (in meV) 
as a function of the wave number q (in units of 10 5 cm -1 ). 
Different line styles refer to different values of the 2DEG 
carrier density rib, ranging from rib = 1 X 10 11 cm -2 to 
rib = 1 x 10 12 cm -2 . The Dirac-fermion density is fixed at 
the value nt = 10 12 cm -2 , while the inter-layer distance is 
d = 30 nm. With reference to Fig. [I] the dielectric constants 
have been fixed to ei = 1 and €2 = 13. Note that A(q) is 
positive, implying a blue shift of the optical plasmon mode of 
the hybrid system with respect to the usual plasmon mode of 
an isolated 2DEG. 



intra-band electron-hole continuum^. Acoustic plasmon 
dispersions lying above the dashed line are not damped 
(within RPA) since they cannot decay by exciting intra- 
band electron- hole pairs in the graphene sheet. When the 
dispersion hits the dashed line, Landau damping occurs 
and the acoustic plasmon pole acquires a finite lifetime. 
Note that when the 2DEG density is low (empty circles 
in Fig. [4]) the acoustic-plasmon dispersion stays out the 
Dirac-fermion intra-band particle-hole continuum for a 
very small range of wave numbers q. It will be thus dif- 
ficult to observe acoustic plasmons in MDSHSs in which 
the 2DEG carrier density is low. 

The acoustic-plasmon group velocity c s = 
lim 9 ^o "^ac (<?) / 9 can be calculated analytically fol- 
lowing the procedure explained in Ref. [2H After 
some tedious but straightforward algebra we find the 
following equation for x — c s /«d, the ratio between 
the acoustic-plasmon group velocity c s and the Dirac 
velocity v^: 



2N t a cc £,Tf(x)g(x) - e 2 Tg(x)\/x 
-2e 2 N t a cc f{x)y/x 2 N h r 2 - Aa. 
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2A b /(fc F , t a B ), f(x) 



x — yx' 2 — 1, 
Note that the 



where £ = e?fcp. t , T = 
and g{x) = x\/N^r s 

solution of Eq. ^ depends only on e 2 , while, as we have 
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FIG. 3: (Color online) Same as in Fig. [2] but for different 
values of the Dirac-fermion density n t and for a fixed value 
of the 2DEG carrier density rib = 2 x 10 11 cm -2 . Notice the 
non-monotonic behavior of A(g): for each q, A(q) reaches its 
minimum when the Dirac-fermion density matches the 2DEG 
density (n t = rib). 



seen above, Lu op (q — > 0) depends on the average e. 

An undamped acoustic plasmon emerges for x > x C rit, 
where the threshold a; cr i t can be directly found from 
Eq. Q. When the Dirac velocity «d is larger [smaller] 
than the Fermi velocity wp,b = ^fcF,b/"ib in the 2DEG, 
x crit = 1 [xcrit = v F ^/vu = 2a cc / {r S ^/N h )]. The ex- 
istence of a solution of Eq. Q thus depends on three 
parameters, namely d, n t , and n b (when the dielectric 
constants e\ and e 2 are fixed). For example, given nt 
and rib, the acoustic plasmon emerges out of the contin- 
uum for d > d( cnt ) with 
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if v D > v F _ h 



otherwise 



(8) 

We remind the reader that the functions f(x) and g(x) 
have been defined right after Eq. ([7]). Similarly, given n t 
and d, the acoustic plasmon emerges out of the contin- 
uum for n > n( CIIt \ where 



(crit) 



e 2 N b (e 2 + 2£T) 

iV b (e 2 4- 2tjN t a cc ) 2 
47re 2 a|a 2 c ( e 2 4- 4£N t a cc ) 



if u D > VF.b 



otherwise 



(9) 

Note that the critical values of d and rib do not depend 
on nt when vr> > «F,b- They indeed depend only on the 
products r/cp,t and £T, which are independent of fcp,t- 
The theoretical prediction based on the solution of 



FIG. 4: (Color online) The acoustic plasmon mode of a mass- 
less Dirac/Schroedinger hybrid system. The acoustic plasmon 
dispersion hwn C (q) (in meV) is plotted as a function of the 
wave number q (in units of 10 5 cm -1 ). Empty symbols label 
data for two different values of the 2DEG carrier density rib. 
In this figure we have fixed n t = 1 x 10 11 cm -2 and d — 60 nm. 



The other parameters are identical to those used in Figs. |2|3| 
In the inset we illustrate the comparison between numerical 
results for the acoustic plasmon dispersion (empty symbols) 
and the analytical result u; a c(g — > 0) = c s q (solid lines) with 
the velocity c s obtained from the solution of Eq. (m). The 
dashed line marks the upper boundary of the Dirac fermion 
intra-band single particle continuum. 



Eq. Q is compared with the numerical results in the 
inset to Fig. [4] 



IV. COULOMB DRAG RESISTIVITY 

In a drag experiment!^ a constant current is imposed 
on one layer (the "active" or "drive" layer) . If no current 
is allowed to flow in the other one (the "passive" layer), 
an electric field develops, whose associated force cancels 
the frictional drag force exerted by the electrons in the 
active layer on the electrons in the passive one. The drag 
resistivity po , is defined as the ratio of the induced volt- 
age in the passive layer to the applied current in the drive 
layer and reads p D = -OD/(ftOt - °d) — _(J D / '(°t0Tt>) 
in terms of the intra-layer conductivities <r t (°b) of the 
top (bottom) layer and the drag conductivity <td • Notice 
that in writing the last equality of the previous equation 
we have assumed an <C at i °V 

Within the Kubo formalism the drag conductivity at 
the second order in the screened inter-layer interaction 
U th {q,Lo) readfP(/i= 1) 
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where (3 = (/cbT) -1 and Ti(q,uj) is the so-called "non- 
linear susceptibility" of the i-th layer. The dynamically- 
screened inter-layer interaction that appears in Eq. (fTol) 



is U t b(q,Lj) = Vtb(g)/e(<7, w), where s(q,uj) is the RPA 
dielectric function dehncd in Eq. 

Here we focus only on the low-temperature limit, i.e. 
k B T -c min(eF,t, £F,b)- In this regime we can do the 
following approximations^: (i) use the low-temperature 
expressions for the intra-layer conductivities <7j, (ii) sub- 
stitute C/tb(g,w) in Eq. (10) with the statically-screened 
inter-layer interaction Uth(q, w = 0), and (hi) expand the 
non- linear susceptibilities Ti(q,uj — > 0) up to the lowest 
order in u). 

Following Ref. |2H we find the following results for the 
top layer (graphene sheet): 



lim cr t = — A r tTt(/cF,t)£F,t 



(11) 



and 
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and the following results for the bottom layer (2DEG): 
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where up.b is the 2DEG Fermi velocity. Note the sin- 
gularity at q — 2/cF,b in the non-linear susceptibility in 
Eq. (14). On the contrary, T t (q,uj —> 0) vanishes at 
q = 2fc F ,t, as a consequence of the absence of backscat- 
tering for massless Dirac fermionJSl 



Using Eqs. ( 11 )-( 14 ) we find that the drag resistivity 



in the low-temperature limit is given by 

i i (fc B r) 2 r 
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\ dq\U th (q,0)\ 2 
(15) 



where q max = min(2fcF )t , 2fc F ,b) an d F{x, y) — 

y/(l-X*)/(l- V 2 ). 

It is convenient at this stage to introduce a dimension- 
less expansion parameter q = d^kp^kF.b- The MDSHS 
is weakly (strongly) coupled if 97 ^> 1 (77 1). Similar 



reasoning to the one employed in Ref. |2H yields the fol- 
lowing result for pn in the weak-coupling limit (restoring 
Planck's constant for clarity): 

h <(3)e 2 (k B T) 2 

hm hm pn = 
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(16) 



We emphasize that the functional dependence of prj on 
the carrier densities n$ and on the inter-layer distance d 
in this limit is the same as in the case of two 2DEGsP 
and two graphene sheetiP^S'. 

In the strong-coupling limit, we find 



lim lim p^ = — 



h 7r 2 v D (A: B T) 2 
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(17) 



where x max = min(2 Ak /A;F l t/fcF.b, 2yJk F ^/k-p^) and we 
have introduced the Thomas-Fermi screening wave num- 
bers g T F,t = N t a cc k Ftt and <j T F,b = N?/ 2 r s k F ^/2. 

The quadrature on the r.h.s. of Eq. ( |17[ ) can be carried 
out analytically in the special case k Ft = ^F.b- Due to 
the different degeneracies in two layers, this condition 
implies a density imbalance between the two layers: rib = 
N h n t /N t = 



lim lim od 

il->OT-sO 



In this case, Eq. (17) yields 

h 7T vo (k B T) 2 
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N t a cc + Nl' 2 r s /2 
21 ) 



2e + N t a cc + N* /2 r s /2 



(18) 



Eqs. (|15|)-(18) are the most important results of this Sec- 



tion^- 1 . As in the case of drag between two graphene 
sheets^!, Eq. ( 18 1 does not depend on the inter-layer dis- 
tance d. In the limit n — > it yields a dependence of prj 
on carrier density of the form pp oc n, -1 , analogously to 
what found by Carrega et alW* for two graphene sheets 
with equal density [note that in the limit n — > one has 
to find the asymptotic behavior of the expression in curly 
brackets in Eq. ( |18[ ) for r s — > 00]. Finally, we emphasize 
that Eq. (18) does not contain the well-known^ hi(T) 



factor which appears in the drag resistivity at strong cou- 
pling in the case of two 2DEGs in semiconductor double 
quantum wells. This factor, which stems from the con- 
tribution to drag from momenta q of the order of 2k F , is 
completely suppressed in the present case by the absence 
of backscattering for massless Dirac fermions^ [compare 
Eq. (pi) with Eq. ((Til)]. 
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V. SUMMARY AND DISCUSSION 

In this work we have proposed that massless 
Dirac/Schroedinger hybrid double layers be used to probe 
Dirac plasmons in a graphene sheet by employing res- 
onant inelastic light scattering^. A natural resonant 
condition in this system is offered by the band gap 
of the semiconductor (e.g. GaAs) hosting an ordinary 
parabolic-band 2D electron gas, thus bypassing the non- 
trivial issue of the absence of a gap in an isolated 
graphene sheet. We have demonstrated that informa- 
tion on Dirac plasmons can be extracted in a rather di- 
rect manner from the measurement of the optical plas- 
mon mode of the hybrid double layer. The latter sup- 
ports also a soft mode which disperses linearly as a 
function of wave number and whose observation requires 
a sufficiently-high electron density in the semiconduc- 
tor quantum well. Finally, we have calculated the low- 
temperature Coulomb-drag resistivity in the Boltzmann- 
transport limit, showing that in the strong-coupling limit 
it displays a dependence on carrier densities that is not 
shared by conventional all-semiconductor double quan- 
tum welLs^. 

The system depicted in Fig. [T] is amenable to exper- 
imental investigations and paves the way for the study 
of many other intriguing phenomena. We expect, for ex- 
ample, that massless Dirac/Schroedinger hybrid double 
layers will display interesting correlated states (and as- 



sociated transport anomalies) induced by inter-layer in- 
teractions deep in the quantum Hall regime. Another 
very appealing subject of investigation could be "hybrid 
exciton condensates". Exciton condensates are elusive 
many-particle systems in which electron-hole pairs con- 
dense below a certain critical temperature in a coher- 
ent superfluicP^. Double-layer structures are extremely 
useful to spatially separate electrons and holes in two 
different layers, thereby suppressing electron-hole recom- 
bination^. Realizing all-semiconductor-based electron- 
hole double layers is a rather difficult taslP^l. On the 
contrary, holes can be trivially induced in the graphene 
sheet depicted in Fig. [I] by gating techniques. Massless 
Dirac/Schroedinger hybrid double layers thus offer un- 
precedented opportunities to realize, probe, and manip- 
ulate novel electron-hole quantum liquids and, hopefully, 
hybrid exciton condensates at sufficiently low tempera- 
tures. 
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